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ABSTRACT 

We present ROSATPSPG observations of the twin-jet radio galaxy 3C 449. The soft X- 
ray emission from this object is dominated by an extended halo with a scale comparable 
to that of the radio source. The asymmetry of the X-ray emission is reflected in that of 
the radio lobes, providing evidence that the behaviour of the jets is strongly influenced 
by the external medium. A region of reduced X-ray surface brightness coincident with 
the southern radio lobe of 3C 449 suggests that the radio source has displaced thermal 
plasma from the X-ray-emitting halo. However, the minimum pressure in the radio 
lobe is considerably lower than our estimates of the pressure in the external medium. 
We discuss the implications for the dynamics of the radio source. 
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1 INTRODUCTION 

3C 449 (B2 2229+39) is a well-known FRI radio galaxy (Fa- 
naroff & Riley 1974) at a redshift of 0.0171. Its symmetrical 
inner jets have been well studied in the radio (Perley, Willis 
& Scott 1979; Cornwell & Perley 1984). On larger scales, the 
southern jet flares into a lobe, while the northern jet con- 
tinues to be collimated until it fades into the noise (Birkin- 
shaw, Laing & Peacock 1981; Jagers 1987; Andernach et al. 
1992; Leahy, Bridle & Strom 1997). 3C449's host galaxy, 
UGC 12064, has a close bright companion (e.g. Balcells et 
al. 1992) and models have been proposed in which the or- 
bital motion of the host is responsible for the symmetrical 
oscillations in the jets (e.g. Lupton & Gott 1982; Hardee, 
Cooper & Clarke 1994). 

UGC 12064 appears in the optical to lie in a poor cluster 
(Zw 2231.2+3732), and the extended optical halo surround- 
ing it and its companion allows it to be classed as a cD 
galaxy (Wyndham 1966). Miley et al. (1983) used the Em- 
stein IPC to show that it is a weak, extended X-ray source, 
and the rotation measure gradients in the jet are thought to 
be produced by the hot gas responsible for the X-ray emis- 
sion (Cornwell & Perley 1984) . In this paper we present new, 
sensitive observations of 3C449 with the ROSAT PSPC, 
which allow us to map in detail the hot gas halo surrounding 
the source. 

Throughout we use a cosmology in which Ho — 50 km 
s _1 Mpc -1 , go = 0. At the distance of 3C449, 1 arcsec 
corresponds to 0.485 kpc. 



2 OBSERVATIONS AND DATA ANALYSIS 

3C 449 was observed for 10 ks with the ROSAT PSPC be- 
tween 1993 Jan 04 and 1993 Jan 10 as part of a programme 
to study the X-ray emission from a sample of radio galax- 
ies drawn from the B2 radio survey as likely candidates for 
the parent population of BL Lac objects (Worrall & Birkin- 
shaw 1994). We analysed the data with the Post-Reduction 
Off-line Software (PROS). After filtering to remove time in- 
tervals with abnormally high particle background, we were 
left with 9.15 ks of good data. 

An extended source was detected with 1840 ± 85 net 
counts in a 10-arcmin-radius circle about the pointing cen- 
tre, corresponding to a count rate of 0.20 + 0.01 s _1 over the 
PSPC energy range (0.1-2.5 keV). The background count 
rate was estimated from an annulus between 10.8 and 16.7 
arcmin from the pointing centre; source and background re- 
gions were corrected for the effects of vignetting using the 
1-keV parametrisation of the vignetting function. A number 
of point X-ray sources, and one extended source, believed to 
be background sources, were excluded from the source and 
background regions. The closest of these, 8.6 arcmin W of 
the source, is coincident with a foreground star; there is no 
bright optical object coincident with the others, though a 
source 15 arcmin to the SE is coincident with a weak radio 
source. 

Fig. [l] shows the extended X-ray structure of the source. 
Extended object and X-ray background analysis software 
due to S. Snowden (Snowden 1997) was used to model and 
subtract non-cosmic background components and correct for 
exposure and vignetting. The resulting 0.4-2 keV image was 
then adaptively smoothed as described elsewhere (Worrall, 
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Birkinshaw & Cameron 1995; Worrall et al. in prep.). The 
central peak in X-ray brightness corresponds to the position 
of the radio core of 3C 449 to within 9 arcsec, consistent 
within the absolute position errors of ROSAT. The central 
regions around the peak are clumpy and elongated in posi- 
tion angle ~ 45° on scales of up to 10 arcmin (300 kpc). The 
reason for this elongation is not clear; there are no bright 
associated galaxies coincident with it on sky-survey plates. 
The X-ray emission on larger scales is asymmetrical, being 
significantly more extended to the south and west than it is 
to the north and east. 

Spectral fits to the data were carried out in the energy 
range 0.2-1.9 keV (SASS channels 6-29 inclusive), for con- 
sistency with earlier work and because this is the energy 
band over which the PSPC point-response function can be 
modelled accurately. A thermal model is a better fit to the 
data than a power law, as expected from the extended na- 
ture of the X-ray emission. The best-fitting single thermal 
model at the redshift of the host galaxy has abundances 100 
per cent solar, kT = 1.211°;^ keV, galactic N H = 5±g X 10 20 
cm -2 , with errors being la for two interesting parameters. 
This value of Nh is consistent with the value interpolated 
from the results of Stark et al. (1992), 1.1 x 10 21 cm" 2 , 
at the la level; fixing galactic Nh at this value and allow- 
ing intrinsic Nh to vary gives a best-fit model with 50 per 
cent solar abundances and kT — 1.14^g Qg keV, with neg- 
ligible intrinsic absorption. More complicated models, such 
as a two-component thermal model or a thermal model plus 
power law, do not improve the fit to the data; there is also 
no significant difference in the temperatures fitted to the 
inner and outer regions of the source consistent with lit- 
tle variation of the gas temperature across the source. The 
count rate and spectrum imply a 0.2-1.9 keV luminosity of 
4.2 x 10 35 W; the source is thus comparable to the brightest 
sources studied by Worrall & Birkinshaw (1994). 

Although the source is not azimuthally symmetrical, 
radial profile fitting allows its size to be roughly charac- 
terised. We therefore extracted a background-subtracted ra- 
dial profile as described by Worrall & Birkinshaw (1994), 
and determined the best-fitting model consisting of a (3- 
model (Sarazin 1986) and an unresolved component. The 
best-fitting models have low /3; with /3 = 0.35 the best-fit 
model had a core radius of 35 arcsec, with a central gas den- 
sity 4.6 x 10 3 m -3 . The radial profile and best-fit model are 
plotted in Fig. |^. The unresolved component contained ~ 2 
per cent of the counts; any nuclear component is therefore 
weak in comparison to the extended emission, as suggested 
by the spectral-model fits. The particular value of (3 chosen 
does not strongly affect the results discussed in the remain- 
der of the paper. Models with f3 < 2/3 formally require a 
cutoff radius to prevent the gas mass becoming infinite, but, 
as noted by Birkinshaw & Worrall (1993), the effect on the 
model emission measure distribution is small over the region 
in which the X-ray-emitting gas is detected. 

3 RADIO- X-RAY COMPARISON 

Fig. | shows the 608-MHz WSRT map of Leahy et al. (1997) 
overlaid on a grey-scale representation of the image of Fig. 
|l[ The southern lobe appears to be embedded in X-ray emit- 
ting material, while the northern plume is clear of it after 



the bend at 5 arcmin (150 kpc) from the core. This may be 
interpreted as evidence for interaction between the outflow 
from the radio source and the ambient medium; the differ- 
ences between the southern lobe and northern plume are 
then attributed to the different external environments, with 
the higher inferred densities to the south being sufficient to 
confine the outflow and form a lobe. 

Further evidence suggesting that the radio emission is 
interacting with the X-ray-emitting medium (rather than 
just being projected on to it) comes from an apparent local 
surface-brightness minimum in the X-ray emission coinci- 
dent with the southern lobe; this suggests that the pres- 
sure from the relativistic particles in the lobe has displaced 
the thermal plasma, which now seems to be forming a rim 
around it. We emphasise, however, that the individual fea- 
tures forming the southern edge of the rim are detected, 
on maps smoothed with a Gaussian with a = 1 arcmin, 
at the 3-5(7 level only; further observations are needed to 
confirm the nature of these features. (We discuss the signif- 
icance of features in smoothed X-ray images in Appendix 
A.) Fig. ^| shows the edge-brightened nature of this rim; 
the background-corrected surface brightness drops by more 
than a factor of 2 between edge and minimum, although 
the cut of Fig. ^| passes through an off-centre blob of X- 
ray emission. When mean count rates are measured from 
the vignetting-corrected 0.1-2.5 keV data in the deficit re- 
gion (defined using the radio contours as a guide) and the 
'rim' (defined as a region extending 2-3 arcmin out from the 
deficit region, but excluding the bright central regions of the 
source), the mean count rate (source and background) in the 
rim [(4.7 ±0.1) x 10" 7 counts arcsec s 1 exceeds that in 
the deficit region [(3.8 ± 0.2) x 10" 7 counts arcsec s 
at the 3a level, so we regard the deficit of X-ray emission 
as significant. Similar X-ray deficits associated with radio 
lobes have been seen in the cluster-centre galaxy NGC 1275 
by Bohringer et al. (1993) and in the powerful FRII radio 
galaxy Cygnus A (Carilli, Perley & Harris 1994). We cannot 
rule out the possibility that the radio lobe coincides with 
the X-ray minimum purely by chance, but the coincidence 
in size between the radio lobe and the spacing of the 'rim' 
features is certainly suggestive. 

The inner and outer regions of 3C 449 have different ro- 
tation measures (RM); Cornwell & Perley (1984) measure 
-210 rad m -2 near the source centre while Andernach et al. 
(1992) find RM ~ -170 rad m -2 in the outer lobes. The bulk 
of the RM is probably contributed by our Galaxy; nearby 
sources have average RM ~ —210 rad m -2 , though the situa- 
tion is complicated by a foreground HI1 region which crosses 
3C449 (Andernach et al.). If the difference between the cen- 
tral and outer regions of 3C 449 is due to the cluster gas, 
then it can be produced with uniform magnetic fields in the 
X-ray emitting gas of roughly 10 pT. 



4 DISCUSSION: DYNAMICS OF THE RADIO 
SOURCE 

If the radio-emitting material has displaced thermal plasma 
then we expect pressures to be comparable in the southern 
lobe and surrounding gas. Figure [s] shows the pressure pre- 
dicted by the best-fitting /3-model (with f3 — 0.35, core ra- 
dius 35 arcsec and kT = 1.14 keV) as a function of distance 
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from the source centre, together with minimum pressures at 
various points in the radio source. It can be seen that if the 
radio source is in the plane of the sky the minimum pressures 
are lower than the inferred hot gas pressure at all points in 
the source, with the discrepancy being greatest (a factor of 
20) at the far ends of the object. Similar results are obtained 
by using other models for the ambient gas: for example, if 
the southern region of the X-ray source is modelled as a uni- 
form spherical shell with inner radius 4.5 arcmin and outer 
radius 6.2 arcmin (Fig. ^), neglecting the off-centre blob of 
emission which may be a foreground or background feature, 
the thermal pressure is 1.7 x 10 -13 Pa, slightly higher than 
the value predicted by the /3-model. 

Underpressuring of the radio lobes with respect to the 
X-ray-emitting plasma has been observed in other low-power 
sources (Morganti et al. 1988; Killeen, Bicknell & Ekers 
1988; Feretti, Perola & Fanti 1992; Bohringer et al. 1993; 
Worrall et al. 1995) . We can attempt to resolve the discrep- 
ancy by considering projection effects on the radio source, 
by considering dynamical pressure balance, by allowing rel- 
ativistic protons to make a significant contribution to the 
energy density, by considering situations where the energy 
densities are far from equipartition, or by allowing inter- 
nal thermal material. These possibilities will be discussed in 
turn. 

The assumption that the source is in the plane of the 
sky may be incorrect. If the source makes an angle to the 
line of sight 8 both the measured distance from the centre 
and the length used in minimum-energy calculation are un- 
derestimated by a factor 1/ sin#; the minimum pressure in 
a given component drops as a result by a factor of roughly 
(sin^) 4 / 7 . Because of this, small angles to the line of sight 

8°) are necessary to achieve pressure balance even in 
the inner regions of the source where the pressures are most 
nearly matched. This is implausible both because of the ob- 
served symmetry of the jets (the degree of symmetry seen in 
high-resolution VLA radio maps would require velocities in 
the jets at 5 kpc from the core of ^ 0.07c, compared with the 
estimated velocities of <; 0.4c at that distance from the core 
in similar objects: Hardcastle et al. 1996, 1997; Laing 1996) 
and because of the implausibly large true linear sizes (<; 5 
Mpc) they would imply for the radio source. This is unlikely 
to provide a solution to the problem, therefore. If we use the 
velocities inferred from observations of other sources, the an- 
gle to the line of sight 9 is predicted from the jet sidedness 
to be ~ 80° . 

Pressure balance could in principle be dynamical if the 
lobe were expanding; if we postulate a spherically symmetric 
wind, driven by the jet, impinging on the external medium, 
then the problem is similar to that of the hot spots of clas- 
sical double radio sources, and the non-relativistic equation 
of pressure balance in the frame of the edge of the lobe is 

Pint + Pint(Vwind — U cxp ) 2 = P cx t + PextVexp 

where u cxp is the outward velocity of the lobe's edge. This 
can be rewritten in terms of the internal sound speed ai nt 
and Mach number M of the wind: 

A >2 2 / Pext ,\ 0?nt , o A ,< if P oxt 1 I 2 

M a int = — 1 — + 2Ma int v cxp + 1 v oxp 

\ -Tint / I y pint J 

(where V is the internal adiabatic index) from which it 
can be seen that no subsonic solution exists if p cx t > pint 



and P ex t S> -Pint- For a purely relativistic internal plasma 
aint ~ c/Vo, so a relativistic wind would be required for 
pressure balance, and we would expect to see effects of beam- 
ing on the surface brightness of the lobe. If cold denser mate- 
rial were mixed into the plasma (e.g. entrained material, dis- 
cussed below) then ajnt would drop and Pi nt and pi nt would 
rise, making lower wind speeds possible. The energetics of 
this process, though demanding, are not completely incon- 
sistent with the power that could be supplied by the jet. It 
is harder to apply this explanation to the underpressuring 
in the inner regions of the source. 

The minimum pressure of the radio plasma is a weak 
function of the ratio k between the number of radiating 
and non-radiating relativistic particles in the plasma, going 
roughly as (1 + re) 4 / 7 . A large contribution from relativistic 
protons (re ~ 200) is necessary before the internal and ex- 
ternal pressures in the outer regions of the source become 
similar. A plasma filling factor <j> of less than 1/200 has the 
same effect. Neither of these situations can be ruled out. Fil- 
amentary structures seen in radio maps of the best-studied 
objects at high resolution may indicate filling factors of this 
order when projection effects are taken into account. In the 
inner jets, where the pressure discrepancy is less, a more 
moderate contribution from protons (re ~ 20) or a larger 
filling factor (<j> ~ 1/20) is sufficient to achieve pressure bal- 
ance. 

The electron and magnetic field energy densities may 
be far from equipartition. The equipartition B-field in the 
southern lobe, treated as a uniform sphere with re = 0, is ap- 
proximately 0.14 nT; to achieve the required pressures, the 
B-field must be either 0.9 nT (a factor of 6 above equipar- 
tition) or 30 pT (a factor of 5 below equipartition). In the 
latter case, the flux in X-rays from inverse-Compton scatter- 
ing of the cosmic background radiation (calculated using the 
code of Hardcastle, Birkinshaw & Worrall 1998) would be 
significant and detectable with our observations, so this pos- 
sibility can be neglected and only a high P-field is consistent 
with observation. However, Feigelson et al. (1995) report on 
X-ray observations of Fornax A which are consistent with a 
lobe magnetic field slightly weaker than equipartition (with 
k = 0, 4> = 1) in that object. Again, in the inner jets, the 
problem is less severe; the P-field must be about a factor 4 
above equipartition, or a factor 3 below it. 

The jets will entrain material as they pass through the 
galaxy (e.g. Bicknell 1994 and references therein), and the 
entrained gas will provide a contribution to the internal pres- 
sure. The fact that the discrepancy between internal and 
external pressures increases with distance from the galaxy 
is consistent with this model, since we would expect the 
amount of entrained material to increase in this way. But the 
fact that there is a deficit in X-ray emission in the S lobe sug- 
gests that there is little hot (kT ~ 1 keV) gas in this region, 
and so that there is little contribution to the internal pres- 
sure from hot gas. Cold gas would not be observed in X-rays 
but would contribute less to Pint- The observed depolariza- 
tion between A = 49 and 21 cm, (DP~ 0.75; Jagers 1987) in 
the centre of the lobe suggests that the internal thermal elec- 
tron density n e ^ 60 m -3 , using the simple models of Cioffi 
& Jones (1980), which is less than the inferred external par- 
ticle density at this radius (n ~ 300 m -3 ). Depolarization, 
at least in the inner regions of the source, can be attributed 
entirely to rotation measure gradients across it (Cornwell & 
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Perley 1984). However, field reversals can hide large quanti- 
ties of thermal material (Laing 1984), so it is possible that 
thermal material at temperatures significantly different from 
1 keV contributes to the internal pressure; the fact that the 
608-MHz polarization is much higher at the edges of the 
lobe than at the centre is certainly suggestive of internal 
depolarization. 

If the radio source is genuinely underpressured with 
respect to the hot gas, the association between the radio 
lobe and the X-ray void must be considered coincidental, 
in the sense that the lobe has simply filled a pre-existing 
region of low gas density. This is not completely out of the 
question, since the collapse time of such a region is com- 
parable to D/cx, where D is the size of the region and 
cx(~ 400 km s _1 ) is the sound speed in the X-ray gas; this 
time is a few x 10 8 years, comparable to the lifetime of the 
radio source. The hypothesis is consistent with the generally 
unrelaxed appearance of the X-ray emission. The problem 
again is to account for the underpressuring in the inner re- 
gions, where the collapse times are much shorter. 



5 CONCLUSION 

We have presented evidence that the large-scale structure 
of 3C 449 is being determined by the distribution of the hot 
intergalactic plasma. 

The radio lobes are significantly underpressured if the 
radio-emitting plasma is electron/positron with filling fac- 
tor unity and the magnetic field and particles are in energy 
equipartition. It is then hard to understand how the south- 
ern lobe can be associated with a deficit of X-ray emission. 
This problem, which seems to be common in low-power ra- 
dio sources, may be resolved in a number of ways; the most 
plausible seem to be a dominant contribution to the parti- 
cle energy density from relativistic or thermal (entrained) 
protons or a plasma filling factor much less than unity. 
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APPENDIX A: UNCERTAINTIES IN 
SMOOTHED X-RAY IMAGES 

It is often necessary to convolve X-ray images with one 
or more smoothing Gaussians in order to analyse extended 
structure. The disadvantage of this procedure is that it ren- 
ders a calculation of the significance of detected structure 
or of the uncertainties on measured quantities more diffi- 
cult. Here we outline a means of overcoming the problem. 

When a two-dimensional array of Poisson noise, with 
mean /i counts per pixel, is convolved with a Gaussian 
(/(r) oc exp(r 2 /2<r 2 )) it is tempting, but incorrect, to treat 
the Gaussian as a single bin of size 27ra 2 and apply Pois- 
son statistics (with a mean of 27rcr 2 /i) to calculate the sig- 
nificance of the results. This approach badly overestimates 
the spread of the cumulative probability distribution (CPD), 
and so overestimates the errors, because it treats the value of 
each pixel in the convolved array as the uniformly weighted 
sum of a small, finite number of single pixels containing 
Poisson-distributed values, rather than as the Gaussian- 
weighted sum of a large (formally infinite) number of such 
pixels. 

A more rigorous treatment produces a better answer. 
Let the number of counts in a pixel in a smoothed image, 
Cj, be the weighted sum of the contributions from the pix- 
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els in the original image, a; Cj = ^WiCi, where Wi are the 
weights. If the expected number of counts per pixel in the 
unsmoothed image was fi, then (cj) = /i^u>i; clearly we 
require J^iUi = 1. The variance of Cj, is given by Var(cj) = 
Wj Var(cj); for Poisson-distributed d, Var(cj) = fi^w 2 . 
If the smoothing function is a 2-D Gaussian adequately sam- 
pled by the pixels, then we can write Wi = Cexp(— r/2a 2 ), 
where the condition J Wi2-Krdr = 1 requires that C = 
1/27TO- 2 ; then 

I* oc 

Var(cj) = fiC 2 / exp(— r/cr 2 )27rrdr 



4ira 2 



that the standard deviation of the distribution is 



\J [i/Ana 2 . This result differs by a factor of only y/2 from 
that derived from the naive approach discussed above. In the 
limit that <r 2 /i 3> 1, this calculation therefore allows an ac- 
curate calculation of the significance of a detection; the Cen- 
tral Limit Theorem guarantees that the cumulative proba- 
bility distribution (CPD) approaches normality, so that the 
statistics can be treated as Gaussian, and the large a means 
that the assumption of a well-sampled convolving function 
is valid. However, when a 2 )j, ^ 1, these assumptions are not 
applicable and application of the result systematically un- 
derestimates the true errors. 

Because of this problem, the most general approach, 
and the one we have adopted in this paper, is to use numer- 
ical (Monte Carlo) methods to determine the correct CPD. 
The simplest numerical route to the CPD involves gener- 
ating a large number of arrays of pixels containing Poisson 
noise with the correct mean number of counts per pixel fi, 
convolving with the Gaussian of interest, and generating an 
approximate CPD from the results. This can then be used to 
find the number of counts per convolved pixel corresponding 
to the 1, 2 and 3-<r probability levels of a normal distribu- 
tion. Results from this procedure are plotted in Fig. Al. As 



expected, the CPD converges on a normal distribution as 
the Gaussian size a tends to infinity and/or for large fi; con- 
sequently, the standard deviation of the smoothed dataset is 
a good estimator of significance level. However, in the more 
usual case in X-ray astronomy where the Gaussian has a a 
of a few pixels and \i <C 1 the noise is far from normally dis- 
tributed, and the standard deviation of the convolved array 
is not a good estimator of the 68 per cent confidence ('lcr') 
level; Monte Carlo methods then seem to be required. The 
approach can easily be extended to calculating the CPD of a 
binned array of smoothed Poisson noise, and automatically 
takes account of the fact that nearby pixels are not inde- 
pendent after smoothing. This method was used to find the 
error bars in Fig. ^. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science P/IpjX style file. 
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Figure 1. Background-subtracted X-ray contours of the extended emission around 3C449 in the PSPC energy band 0.4-2 keV (SASS 
channels 13 to 30 inclusive), adaptivcly smoothed as described in the text. Contour levels (above the background) are 0.15 X 10~ 3 X 
(1, 2, 3, ... , 10, 15, 20, . . . , 50) counts arcmin 2 s -1 . (The corresponding negative contours enclose no data and so do not appear.) The 
cosmic-background contribution was estimated to be 0.27 X 10~ 3 counts arcmin 2 s _1 . The large cross marks the optical position of 
3C449 and the smaller crosses those of two nearby galaxies. 
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Figure 2. Background-subtracted radial profile of the X-ray emission around 3C449 between 0.2 and 1.9 keV. The dotted lines show 
the two components of the best-fit model and the dot-dashed line shows the total model for comparison with the data. \ 2 = 7.7 with 11 
degrees of freedom. 
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Figure 4. East- west slice through the X-ray emission from 3C 449 
at the widest point of the southern 'ring' of X-ray emission, cen- 
tred on RA 22 28 54.4 DEC +38 58 54. West is to the left. 0.1-2.5 
keV X-ray data (all channels) binned into 15-arcsec pixels and 
smoothed with a Gaussian with a = 1 arcmin; each data point 
represents the average over a bin with width 0.5 arcmin (mea- 
sured east-west) and length 2.5 arcmin (measured north-south). 
Superposed is the profile expected from a spherical shell with 
inner radius 4.5 arcmin and outer radius 6.2 arcmin, arbitrarily 
normalised. The dotted line shows the background level measured 
from the background annulus described in section ^; bars cross- 
ing it mark the approximate size of the southern radio lobe. Er- 
ror bars, spaced at one per 4 data points, are ±1<t (68 per cent 
confidence level) for the given background noise level, convolving 
Gaussian and binning scheme, derived from simulation (see Ap- 
pendix A). The off-centre excess in the data is interpreted in the 
text as foreground or background emission. 
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Figure 5. Pressure from the hot gas as a function of distance 
from the centre, calculated from the best-fit /3-model with core 
radius 35 arcsec and /3 = 0.35. Superposed are measurements of 
minimum pressure of various regions in the radio source. Squares 
denote measurements made north of the core, circles those made 
to the south. 



© 0000 RAS, MNRAS 000, 000-000 



10 Hardcastle, Worrall and Birkinshaw 



CD 




Gaussian size (pixels) Gaussian size (pixels) Gaussian size (pixels) 

Figure Al. Effects of different convolving Gaussians on the '3cr' level of an array of initially Poisson-distributcd noise. The level of a 
contour that would exclude all but 0.135 per cent of the background noise, in counts pixel -1 , is plotted (solid line) as a function of the 
size in pixels of the convolving Gaussian for three different choices of the original mean number of counts per pixel (fi). On the left, 
fj, = 0.01 counts pixel -1 . In the centre, fi = 0.1 counts pixel -1 . On the right, fi = 1.0 counts pixel -1 . Plotted for comparison (dashed 
line) is the level that would be correct if the statistics were Gaussian (mean plus 3 times the standard deviation) . It will be seen that the 
two are increasingly similar for large convolving Gaussian and large /i, as expected. Dotted lines in the left-hand figure show the effects 
of assuming that convolving with a Gaussian is equivalent to binning with a bin size equal to the effective area of the Gaussian; they 
represent the upper and lower bounds from Poisson statistics on the number of counts corresponding to the 99.865 per cent confidence 
level. As discussed in the text, this approach systematically overestimates the errors. 
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